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Collectivity in the light radon nuclei measured directly via Coulomb excitation 
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Background: Shape coexistence in heavy nuclei poses a strong challenge to state-of-the-art nuclear models, 
where several competing shape minima are found close to the ground state. A classic region for investigating 
this phenomenon is in the region around Z = 82 and the neutron mid-shell at A = 104. 

Purpose: Evidence for shape coexistence has been inferred from o-decay measurements, laser spectroscopy 
and in-beam measurements. While the latter allow the pattern of excited states and rotational band struc¬ 
tures to be mapped out, a detailed understanding of shape coexistence can only come from measurements 
of electromagnetic matrix elements. 

Method: Secondary, radioactive ion beams of ^®^Rn and ^®®Rn were studied by means of low-energy 
Coulomb excitation at the REX-ISOLDE facility in CERN. 

Results: The electric-quadrupole {E2) matrix element connecting the ground state and first-excited 
2^ state was extracted for both ^®®Rn and ®®®Rn, corresponding to B(E2-,2^ —^ 2]®) = 291*^® W.u. 
and iStll W.U., respectively. Additionally, E2 matrix elements connecting the 2f state with the 4)® 
and 2^ states were determined in ®®®Rn. No excited O’® states were observed in the current data set, 
possibly due to a limited population of second-order processes at the currently-available beam energies. 

Conclusions: The results are discussed in terms of collectivity and the deformation of both nuclei studied 
is deduced to be weak, as expected from the low-lying level-energy schemes. Comparisons are also made 
to state-of-the-art beyond-mean-field model calculations and the magnitude of the transitional quadrupole 
moments are well reproduced. 
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Shape coexistence in nuclei is a phenomenon whereby 
two or more nucleon configurations, each with a differ¬ 
ent macroscopic shape, exist together at similar energy. 
It has been observed in a number of regions of the nu¬ 
clear chart and, over the last decade and more, extensive 
experimental evidence has been found in support of the 
shape coexistence in the Pb region [1]. The most-striking 
early indications came from isotope-shift measurements 
in mercury (Z = 80), which showed a large discontinu¬ 
ity in the mean-square-charge radii between ^®®Hg and 
is^Hg [2]. This was interpreted as a dramatic change 
in shape using calculations based upon the Strutinsky 
shell-correction method [3]. The ground states of the 
heavier isotopes were calculated to be weakly deformed 
and oblate in nature, but when approaching the neutron 
mid-shell at N = 104, this picture changed to a more- 
strongly deformed prolate shape. These shapes are asso¬ 
ciated with structures based upon two different proton- 
hole excitations across the Z = 82 shell closure, namely 
7r(0p-2/i) and 7r(2p-4/i). Recently, the first direct evidence 
of shape coexistence in the even-mass Hg isotopes came 
from Coulomb-excitation experiments [4, 5] , which quan¬ 
tified the deformation of ground and excited 0+ states for 
the first time in this region. 

At Z = 82, the lead isotopes remain spherical in 
their ground state all the way to mid-shell, as indi¬ 
cated by isotope-shift measurements employing laser 
spectroscopy [6, 7]. In the case of ^®®Pbio 4 , competition 
between three shape minima is observed - oblate, prolate 
and the spherical ground state. This was inferred from a- 
decay measurements of ^®°Po [8], where the three states 
lying lowest in energy were observed to be O'*" states. This 
triple shape coexistence is apparent all around the mid¬ 
shell in the parabolic behaviour of the intruder energy 
levels as a function of mass number [Figure 3 of Ref. 9], 
recently investigated down to ^^“Pbgs [9]. 

The phenomenon persists in nuclei above Z = 82, 
where the polonium isotopes were recently observed to 
have a much earlier and more gradual onset of deforma¬ 
tion than observed in mercury [10], without the unusual 
odd-even staggering [11]. One might consider that mer¬ 
cury {Z = 80) and polonium {Z = 84) are analogues with 
respect to their nucleon configuration; the oblate struc¬ 
ture in the mercury isotopes, driven by 7r(0p-2/i) config¬ 
urations, should manifest itself in polonium in 7r(2p-0/i) 
configurations and similarly for the prolate structure. In¬ 
deed, the same parabolic behaviour of intruder states was 
observed when approaching mid-shell [12] and was inter¬ 
preted to be of 7r(4p-2/i) configuration [13]. Coulomb- 
excitation measurements recently determined multiple 
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FIG. 1. (Color online) Systematics of excited states in even- 
even radon nuclei. The lowest-known even-spin positive- 
parity states are shown, with the first- (red) and second- 
excited (blue) states each connected by a solid line to guide 
the eye). 


low-lying matrix elements for nuclei in the transitional 
region where the onset of deformation is observed [14]. 

The radon isotopes {Z = 86) can be expected to have 
similar proton-hole analogues to the platinum isotopes, 
where spectroscopic information on deformed intruder 
states exists beyond the neutron mid-shell [15]. The 
energy level systematics of the even-spin positive-parity 
states in the light even-mass radon isotopes are shown 
in Figure 1, where one can observe decreasing excita¬ 
tion energy of the 2+ state towards ^®®Rn [16, 17]. A 
corresponding deviation from sphericity at = 116 is 
observed in the mean-square-charge radii [18, 19], earlier 
still than in the Po isotopes [10]. This may indicate that 
there is indeed a region of deformation towards the neu¬ 
tron mid-shell that is unreachable within the current ex¬ 
perimental limitations. A more detailed understanding, 
with complementary experimental probes, of the isotopes 
around this transition region, ^®®“^°^Rn(7V = 112—118), 
would help to determine if this behaviour is in fact due 
to the presence of shape-coexisting intruder states. 

Low-lying excited states in the isotopes around the 
N = 126 shell closure are generally considered to be 
associated with a seniority scheme [20, 21], while lower 
masses are proposed to be candidates for vibrational nu¬ 
clei. The observed equal level spacing i.e. a ratio of the 
4+ to 2+ excitation energy (R 42 ) close to 2, indicates a 
possible vibrational nature and the existence of a second 
2~^ state at a similar energy to the 4^ in ^°^Rn adds fur¬ 
ther weight to this argument. The harmonic quadrupole 
vibrator should lead to a very definite and simple pat- 
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tern of states with a single-phonon state with = 2+, 
a triplet of two-phonon states with = 0+, 2+ and 4+, 
and so on. As far as ^°^Rn and ^°^Rn [22] are concerned, 
several of the expected members of vibrational multiplets 
are missing, although it is not presently clear if this is 
due to an experimental limitation. Their low-lying level 
schemes of interest to this study are shown in Figure 2. 
In particular, no observations of a excited 0"^ states have 
been made in any of these nuclei. Although its existence 
is expected in both a vibrational and intruder picture, 
the energy and B{E2) values connecting 2+ states would 
definitively determine the structure. Detailed investiga¬ 
tions of excited 0+ states into the cadmium isotopes have 
proceeded in a similar vein [23], where the vibrational 
picture was found not to be adequate [24] . Additionally, 
the presence of two near-parallel fi"*" states in both nuclei 
is hard to accommodate in a simple vibrational picture. 

A detailed understanding of shape coexistence, or vi¬ 
brational nuclei, will never be achievable from a sin¬ 
gle class of measurement. A comprehensive picture of 
the underlying physics can only come from extraction 
of electromagnetic matrix elements involving a comple¬ 
mentary set of experimental probes. Transition matrix 
elements may be derived from lifetime measurements, in 
combination with precision branching and mixing ratios, 
from in-beam or decay spectroscopy. Coulomb excita¬ 
tion allows not only the extraction of transition matrix 
elements but also of diagonal matrix elements, includ¬ 
ing their sign. These can be used to further conclude 
on the sign of the spectroscopic quadrupole moment for 
excited states and hence, the type of nuclear deforma¬ 
tion. Multi-step Coulomb excitation needed to investi¬ 
gate low-lying non-yrast states in these nuclei, requires 
the availability of intense accelerated radioactive ISOL 
beams, which have only recently become available at fa¬ 
cilities such as SPIRAL and REX-ISOLDE. A pioneer¬ 
ing example of this technique was in [25] at SPI¬ 


RAL. An intense program of Coulomb-excitation experi¬ 
ments at REX-ISOLDE has been underway to study the 
Z = 82 region. This facility is chosen as it is uniquely 
capable of producing beams of heavy proton-rich nuclei 
from spallation reactions. Eurthermore, key techniques 
such as laser ionisation have been developed to produce 
isobarically-pure secondary beams. Experiments involv¬ 
ing very heavy (A > 200), post-accelerated beams have 
proven successful at REX-ISOLDE in recent years, in¬ 
cluding those employing radon [26]. Studies such as 
these, performed at ISOL facilities around the world, are 
currently pushing the boundaries of nuclear spectroscopy 
on the precision frontier in exotic nuclei [27]. 

In addition to the possibility of measuring electro¬ 
magnetic matrix elements. Coulomb excitation is a well 
adapted technique for locating missing states, especially 
low-lying, non-yrast states that may not otherwise be 
populated in-decay or fusion evaporation experiments. 
Since low-lying 0'*' states are key to the understanding 
of these nuclei, exploring the possibility of populating 
a 0^ state via a two-step Coulomb-excitation process is 
desirable. 


II. EXPERIMENT AND DATA ANALYSIS 

Radioactive beams of ^°^Rn and ^°^Rn were produced 
at the ISOLDE facility in CERN via bombardment of 
a uranium-carbide primary target with 1.4-GeV protons 
from the PS Booster. The target-ion-source coupling in 
this experiment was key to reduce isobaric impurities ex¬ 
pected when working with a noble-gas beam. A plasma 
ion source [28] was utilised and an extraction voltage 
of 30 kV was applied along the transfer line and con¬ 
tinuously cooled by a water flow in order to suppress 
the transport of less volatile elements. At the begin¬ 
ning of the running period, the yield of the two radioac¬ 
tive species were measured using the dedicated ISOLDE 
tape station and found to be 9 x 10® ions//rC (^°^Rn) 
and 2 x 10^ ions//rC for ^®"^Rn. The singly-charged 
ions were accumulated and cooled in an ion trap, REX- 
TRAP [29, 30]. At intervals of 58 ms, the potential bar¬ 
rier was lowered allowing bunches of cooled ions to escape 
into an electron-beam ion source, REX-EBIS [29, 30], 
where the charge state of the ions was increased by charge 
breeding up to 47“''. The ^°^Rn and ^°'^Rn beams were 
then accelerated to 2.9 and 2.845 MeV/u, in the 2008 
and 2010 campaigns, respectively, by the REX linear ac¬ 
celerator [31]. A failure of the 9-gap resonator, the final 
element of the REX-LINAC, in the original 2008 cam¬ 
paign restricted the running period. This meant that a 
significant amount data, for both ^°^Rn and ^°^Rn, was 
taken at the lower beam energy of 2.28 MeV/u. 

The secondary radioactive beams were incident on thin 
metallic foil targets positioned at the centre of the Mini¬ 
ball germanium detector array [32]. The delivered beam 
currents at the target position were estimated to be 
around 3 x 10^ ions/s for ^®^Rn and 2 x 10® ions/s for 
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FIG. 3. (Color online) A two-dimensional spectrum, as a 
function of energy and laboratory angle, of scattered projec¬ 
tiles and recoils in the CD detector in the ^°®Ag(^°^Rn) re¬ 
action, at a beam energy of 2.85 MeV/u and target thickness 
of 1.9 mg/cm^. There is no condition on the detection of a 7 
ray. 


The isobaric purity of the beam was monitored 
through inspection of the 7 -ray spectrum obtained with 
a germanium detector positioned at the beam dump, 
approximately 3 m downstream of the target chamber. 
Aside from transitions due to normal room background, 
the beam-dump spectrum only contained 7 -ray transi¬ 
tions following the /3- and a-decay of the isotopes of in¬ 
terest. However, during the second campaign in 2010, 
the cooling of the transfer line of the ion source failed. 
This higher temperature allowed volatile elements to 
pass, in particular, a significant amount of the stable 
^°^Hg, which caused contamination of the beam. From 
the observation in the Coulomb-excitation spectrum of 
the — 7 ^ Oj*") transition at 439.5 keV, in com- 

bination with the previously measured B{E2\2^ —>■ 0 )*') 
value [33], the integrated beam current associated with 
202 Hg was deduced and represented 10 % of that associ¬ 
ated with ^°^Rn. 

The Miniball array [32] comprises eight triple-cluster 
germanium detectors; each crystal is six-fold segmented, 
leading to a total of 144 discrete detector elements. The 
total efficiency of the array is « 7% for 1.33 MeV 7 
rays. Scattered heavy ions were detected in an 500-^m- 
thick annular silicon double-sided silicon strip detector 
(DSSSD) segmented into four quadrants. This CD detec¬ 
tor, so-called due to its resemblance to a Compact-Disc, 
has 16 annular strips on the front face and 24 sectors on 
the back, and covered the range of laboratory angles from 
9 = 16.2° to 53.3°. Figure 3 illustrates the kinematics 
for scattering of ^°^Rn on ^^^Ag at a centre-of-target en¬ 
ergy of 535 MeV. The reactions are performed in inverse 
kinematics so there are two solutions for the projectile 
case. In addition, for the lowest laboratory angles, there 
is an ambiguity between the scattered projectile and tar¬ 
get ions, such that the first two strips of the CD detector 
cannot be utilised in the analysis. In order to resolve the 
issue of having an ambiguous conversion from laboratory 
angle to centre-of-mass (CoM) scattering angle, crucial 


for the calculation of the Coulomb-excitation cross sec¬ 
tions, a coincidence gate on the recoiling target nuclei is 
applied. Here, the events corresponding to the second 
solution, at the very lowest CoM scattering angles, are 
not detected since the recoils do not have enough energy 
to exit the finite width of the target from the point of 
reaction. Any of those that do (e.g. when the reaction 
occurs at the back of the target) are below the energy 
threshold of the CD detector. Therefore, one can confi¬ 
dently assume that all recoil events are from the higher 
CoM-scattering-angle solution. 

Triggered by the release of EBIS, data is collected from 
all detectors during 800-/xs wide “beam-on” window fol¬ 
lowing by an equally-wide “beam-off” window 4-10 ms 
later. In software, a correlation window of 6 /xs is defined 
around the each 7 -ray event of the “beam-on” window, 
and all particles that fall within this window are asso¬ 
ciated with that 7 -ray. In this way, it is possible for a 
single particle to be correlated to multiple 7 rays, but not 
vice versa. Prompt and random windows are defined by 
taking the time difference between the particle and 7 -ray 
triggers, as shown in Figure 4. The particle multiplicity, 
shown in Figure 5, of each event can now be defined as 
mp-nr, where m(n) is the number of prompt (random) 
particles. In order to subtract the randomly coincident 
background from the spectra, 0p-2r and Op-lr events are 
treated in exactly the same way as their prompt coun¬ 
terparts, but given a weight of —Tp/T^, where is the 
width of the prompt and random time windows, respec¬ 
tively. The 7 -ray spectra of Figures 6 , 7, 8 , and 9 show 
background-subtracted 2 p- 0 r events, where each parti¬ 
cle can be identified as a recoil and projectile coincident 
within a particle-particle time window of 150 ns. Ad¬ 
ditionally, Ip-Or events, where only the recoil is identi¬ 
fied, are included. In these cases, the projectile kine¬ 
matics, i.e. laboratory angle and exit energy, must be 
reconstructed for the purposes of performing an optimal 
Doppler correction. Here, two-body elastic scattering is 
assumed and the energy loss is calculated by integrating 
phenomenological stopping power curves fitted to data 
from SRIM [34]. By definition, the 7 -ray multiplicity is 
one, but 7-7 events can still be built by looking for events 
that have the same prompt particle correlations. In these 
experiments, such 7-7 coincidences didn’t provide any 
additional information (see inset of Figure 9). As shown 
in Figure 5, higher-order particle multiplicities do not ac¬ 
count for a significant fraction of the data and are not 
taken into account in this analysis. Events where both 
m and n are greater than zero, i.e., at least one prompt 
particle and at least one random particle, are also not 
considered due to ambiguity in assigning prompt or ran¬ 
dom status. In the case that this represents a significant 
amount of data, it is possible to assume a prompt nature 
for such events, but the weighting of random events must 
be re-considered to account for this. Usually, the ratio 
of intensities of transitions associated with ,5-decaying 
daughter products of the beam, assumed to be purely 
random in time, is then used. 
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FIG. 4. Time difference between a 7 ray, which is acting 
as a trigger, and all correlated particles. The effects of the 
800-ns downscaling window is clearly visible. Two regions 
are indicated in order to define a particle as in “prompt” 
coincidence, or in “random” coincidence. 



FIG. 6 . (Golor online) 7 -ray de-excitation spectra asso¬ 
ciated with the Goulomb excitation of ^°^Rn on ^°®Ag at 
2.90 MeV/u, Doppler-corrected for projectiles (black) and 
target recoils (red). Only events identified in prompt coin¬ 
cidence with a recoiling target nucleus are shown; random 
events, with respect to the particle -7 coincidence time, have 
been subtracted. Peaks are marked with their energy in keV. 
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FIG. 5. Particle multiplicity for each 7 -ray: the number of 
particles that are defined as prompt or random according to 
Figure 4. 
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FIG. 7. (Golor online) 7 -ray de-excitation spectra asso¬ 
ciated with the Goulomb excitation of ^°'*Rn on ^°®Ag at 
2.90 MeV/u, Doppler-corrected for projectiles (black) and 
target recoils (red). Only events identified in prompt coin¬ 
cidence with a recoiling target nucleus are shown.; random 
events, with respect to the particle -7 coincidence time, have 
been subtracted. Peaks are marked with their energy in keV. 


Due to the inherent difficulties in performing an abso¬ 
lute normalisation to elastically-scattered particles with 
Miniball, caused by an imprecise knowledge of the dead¬ 
time with different coincidence conditions, normalisa¬ 
tion to the excitation of the target is preferred [35]. 
In this case, the ^°^Rn(^°"‘Rn) beam was incident on 
a 4.0(1.9) mg/cm^ target of ^°®Ag, for which the rel¬ 
evant matrix elements are sufficiently well-established 
experimentally. The resulting de-excitation 7 -ray spec¬ 
tra are shown in Figures 6(7). In the excitation pro¬ 
cess on the ^°®Ag target, only the states are popu¬ 
lated in the ^°^’^°^Rn projectiles. This means that the 
B{E2; 2f —>■0]'’) and Qsi^t ) can be determined by utilis¬ 
ing the first-order assumption that matrix elements con¬ 
necting higher-lying states, of which we have no direct 
experimental information, do not contribute. 

Due to the presence of de-excitation 7 -rays from the 


target that are Doppler shifted differently to those from 
the projectile, it is not easy to locate weak 7 -ray transi¬ 
tions in the projectile. Accordingly, data was also taken 
on a 2.0 mg/cm^ target of ^^°Sn, chosen to reduce the 
number and intensity of 7 -ray transitions resulting from 
target excitation. A high-lying first-excited 2+ state at 
1171 keV, with a relatively small B{E2) of 11.4 W.u., 
means that this state is not strongly populated. This re¬ 
duces the complexity of the spectra as well as the back¬ 
ground from Compton-scattered, escaped events, as can 
be seen in Figures 8 and 9. In all of the 7 -ray spec¬ 
tra, the intensity of radon K X-rays is markedly high, 
inconsistent with the expectations of internal conversion 
of E2 transitions. The residual fraction of these X rays 
is associated to K-vacancy creation in atomic collisions 
between the high-Z beam and target [36]. Population of 
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FIG. 8 . (Color online) As in Figure 6 but for the ^^°Sn target. 
The inset shows an expanded portion of the spectrum, with 
a bin width of 4 keV. Peaks are marked with their energy in 
keV. 



FIG. 9. (Color online) As in Figure 7 but for the ^^°Sn target. 
Peaks are marked with their energy in keV. The inset shows 
the background-subtracted 7-7 matrix, gated on the 543-keV 
2f —> 0)*" transition, with a bin width of 4 keV. 


higher-lying states in ^^"‘Rn was inconclusive based on 
the ^^°Sn-target data (see Figure 9), possibly due in part 
to a poor 7 -ray resolution caused by noise on the CD 
detector in this part of the experiment, which affected 
the Doppler correction. It may also be that the popula¬ 
tion of the states was simply below the detection limit 
of the experiment. An upper limit for the observation 
of the 4f —>■ 2'^ transition was determined. In ^°^Rn, 
the 4^ state is a little lower in energy and there exists a 
previously-observed 2^ state at 1029 keV. Both of these 
states are clearly populated in the ^^°Sn-target data (see 
Figure 8 ), albeit with low intensity. It helps that both 
transitions sit at higher energy in the spectra than the 
dominant 2^ 0 )'' transitions, since they are clear of the 

Compton background and can be fitted with a smaller 
uncertainty. The extracted intensities are presented in 
Tables I and II. No additional data is obtained from the 
lower-beam-energy runs and it is not considered in the 
cross-section analysis due to large uncertainties on 7 -ray 
intensities. 

Aside from the known 4^ state in ^°"^Rn, there is the 


TABLE I. Intensities of 7 -ray transitions observed in the cur¬ 
rent Coulomb-excitation experiments of ^°^Rn. Efficiency cor¬ 
rection has been performed, relative to the 2^ —> 0 )^ transi¬ 
tion of the projectile in each experiment. Transition intensi¬ 
ties in the Ag target are also included and can be identihed 
by the odd-spin transitions. 


Beam Energy 

Target 

Transition 

I, 

2.845 MeV/w 


2 +^ 0 + 

990(37) 


( 2.0 mg/cm^) 

2 +^ 2 + 

27(7) 



4+^2+ 

29(6) 



2 + -> 0 + 

< 19(9) 

2.90 MeV/w 

109 Ag 

2 +^ 0 + 

923(40) 


(4.0 mg/cm^) 

3/2/ ^ 1/2/ 

1260(60) 



5/2/ ^ 1/2/ 

1000(50) 


TABLE 11. Intensities of 7 -ray transitions observed in the 
current Coulomb-excitation experiments of ^°'*Rn. Efficiency 
correction has been performed, relative to the 2 )^ —>■ 0 )^' tran¬ 
sition of the projectile in each experiment. Transition intensi¬ 
ties in the Ag target are also included and can be identihed 
by the odd-spin transitions 


Beam Energy 

Target 

Transition 

I, 

2.845 MeV/w 

i"“Sn 

2 +^ 0 + 

6130(200) 

( 2.0 mg/cm^) 

2 + -> 2 + 

< 190(160) 



4+^2+ 

< 240(90) 

2.90 MeV/u 

109 Ag 

2 +^ 0 + 

660(40) 

(1.9 mg/cm^) 

3/2/ ^ 1/2/ 

720(50) 



5/2/ ^ 1/2/ 

700(50) 


potential for the population of an unobserved 2 ^ state. 
Assuming it decays predominantly to the 2^ state as in 
^°^Rn, the 7 -ray de-excitation could form a doublet with 
the 2'^ —>• Oj*" transition at 543 keV. This would place the 
state around twice the energy of the 2^ state, something 
that is expected with a vibrational-like structure. To in¬ 
vestigate this possibility, the background-subtracted 7-7 
matrix for the ^^°Sn-target data was projected with a 
gate between 520 and 570 keV, as shown in the inset of 
Figure 9. A Icr upper limit of a peak [37] was deter¬ 
mined for the region between 520 and 570 keV of 21(19) 
counts. The 7-7 efficiency, was determined at 

311 keV through the 5/2—>• 3/2/ —>•1/2/ cascade in 
^°®Ag and extrapolated using the singles efficiency de¬ 
termined for a ^®^Eu/^^^Ba source combination to give 
e.y.y(543 keV) = 11(3)%. Consequently, the Icr upper 
limit of the number of counts in the singles spectrum is 
190(160) counts, which corresponds to less than 3% of 
I'y(2'^ —>0/") transition. Assuming a similar excitation 
probability for the Ag (Z = 47) target as the Sn {Z = 50) 
target, one can assume that this would not significantly 
affect the determination of the B{E2;2^ —>0/") value, 
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since it is less than the statistical precision of the transi¬ 
tion intensity. 


III. RESULTS 

For the Coulomb-excitation analysis, the GOSIA [38, 39] 
code was utilised in order to calculate excitation proba¬ 
bilities, and consequently de-excitation y-ray intensities, 
for a given set of electromagnetic matrix elements. The 
calculated intensities are then compared to experimen¬ 
tal data, along with additional spectroscopic information, 
such as excited-state lifetimes, E2/M\ mixing ratios and 
y-ray branching ratios. Conversion coefficients used in 
GOSIA were calculated using the Bricc data tables [40]. 
A x^-like, least-squares function is constructed and can 
be minimised with respect to the electromagnetic matrix 
elements as input parameters, along with a set of nor¬ 
malisation constants. For the cases where normalization 
to the ^°®Ag target excitation was used, a special version 
of the code, G0SIA2, is employed. Here, the total is 
calculated for fixed values of the projectile matrix ele¬ 
ments, (0[j''l|i?2jj2[j^) and {2'^\\E2\\2'^), scanning a large- 
scale two-dimensional surface in order to search for the 
best solution at Xmin- The associated Icr uncertainties 
can then be extracted by cutting the surface at Xmin + 1 
and projecting the limits to the relevant axis. These pro¬ 
cedures are described in detail in Ref. [35]. 

In the first step, the level schemes as shown in Fig¬ 
ure 2 are defined in GOSIA2, where the bj*" and 4^ are 
buffer states in ^°^Rn (4]^ in ^°‘^Rn) to prevent an artifi¬ 
cial build-up of population in the highest-energy observed 
states. Both the ^°^Rn and ^°'^Rn data are segmented 
into five different angular ranges, utilising the segmenta¬ 
tion of the CD detector, yielding five independent experi¬ 
ments. This gives a total of five data points in the projec¬ 
tile system; the intensity {I^) of the 2'^ —>■ 0]^ transition 
in each experiment. There are, however, seven parame¬ 
ters; the matrix elements (0]''j|i?2jj2[j^) and (2]'"j|i?2jj2[j^), 
plus five normalisation constants, which can be consid¬ 
ered as a product of the integrated beam current, live¬ 
time fraction and particle-y efficiencies at the given scat¬ 
tering angle. The target system is over-determined with 
ten transition intensities, two in each of the five inde¬ 
pendent experiments (angular ranges), in addition to the 
nine additional spectroscopic data presented in Table III, 
fitted to a total of seven matrix elements and five normal¬ 
isation constants. These five normalisation constants are 
shared between the projectile and target systems and can 
be fitted simultaneously in both data sets. This allows for 
an over-determination of the whole system, which can be 
reduced to a two-parameter system with five data points 
for the projectile. 

Two-dimensional surfaces are plotted in Figures 10 
and 11 for ^°^Rn and ^^‘^Rn, respectively. The kinematics 
of the experimental set-up limited the observed range of 
CoM scattering angles, which in turn limited the sensitiv¬ 
ity to the diagonal matrix element, (2]''j|i?2jj2[j^). A good 


TABLE III. Spectroscopic data related to the low-lying level- 
scheme (1/27, 3/27, 5/27) of ^°®Ag included in the GOSIA 2 
fit. An average was taken of the two possible solutions for 
(5/27IIA 2 || 5 / 27 ). Matrix elements connecting the higher- 
lying states (3/27, 5/27) were determined from previous 
Coulomb-excitation measurements [ 41 , 42 ] and fixed in the 
fit. The 9 / 2 + isomeric state was not included. 

Spectroscopic data for ^°®Ag Reference 


B(£;2;l/27 3/27) 

0.222(19) e^b^ 

[41] 

B(£;2; 1/27 ^5/27) 

0.320(26) e^b^ 

[41] 

(3/27||E2||3/27) 

-1.3l° j eb 

[42] 

(5/27||E2||5/27) 

-0.21 or -0.56 eb 

[43] 

G (5/27 ^3/27) 
7^(5/27^1/27) 

0.069(16) 

[44] 

r(3/27) 

8.5(10) ps 

[45] 

r(5/27) 

47(2) ps 

[46] 

5(3/27 ^ 1/27) 

-0.196(27) 

[41] 

5(5/27 ^ 3/27) 

-0.039(17) 

[41] 
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FIG. 10. (Color online) Two-dimensional total-x^ surface for 
^°^Rn on ^°®Ag at 2.9 MeV/u, extracted from GOSIA. The 
data were segmented in to five angular ranges. 


determination of the spectroscopic quadrupole moment, 
Qs{2,i), which is proportional to {2^^^\\E2\\2'^), requires 
not only statistical precision but data at both high and 
low scattering angles. This in turn achieves a variation 
in sensitivity to subtle higher-order effects. For the cur¬ 
rent data set, a strong overlap in the x^ functions of the 
different experiments leads to an elongation of the la 
confidence region in the (2+||i72||2j^) axis. The strong 
correlation between the two parameters means that the 
determination of (0]''||i72||2[j'") or B{E2;2^ 0]*') is also 

adversely affected, increasing the projected uncertainty. 
Under the assumption of no second-order effect for Qs, 
the uncertainty is equivalent to that of the statistical un¬ 
certainty of the y-ray intensity, but under-estimates the 
true uncertainty by a factor ~ 3.5 in both ^°^’^°‘*Rn. 

Following the extraction of B{E2;2^ —)■ Oj*"), the 
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FIG. 11. (Color online) Two-dimensional total-x^ surface for 
^°^Rn on at 2.9 MeV/u, extracted from GOSIA. The 

data were segmented in to five angular ranges. 


TABLE IV. Transition strengths, B{E2), and spectroscopic 
quadrupole moments, Qs, along with their uncertainties ob¬ 
tained from the two-dimensional analyses and GOSIA min¬ 
imisation. In the case of ^°^Rn, the final values are extracted 
from the full simultaneous analysis of data on both the Ag 
and targets. The uncertainties include correlations to 

all seven matrix elements in the fit. The fit is shown to con¬ 
verge with the two-dimensional x^ analysis and produces con¬ 
sistent uncertainties, proving that the correlations are small. 


■‘Rn 


B{E2-,2+ ^ 

^0+) 

29l® W.u. 

43l]^ W.u. 

B{E2-2+^, 

^2+) 

160+®° W.u. 

- 

B{E2-2+^Q+) 

< 0.4(3) W.u. 

- 

B{E2-4t ^ 

^2+) 

63t]® W.u. 

< 74(30) 

Qs{2+) 


0.9l®;® eb 

-0.4l];° eb 


higher-statistics data for Coulomb excitation on the 
^^°Sn target is analysed in a second step. The ratio 
of transition intensities —>■ “^0^) can 

be considered to be almost directly proportional to the 
B{E2]A^ —>■2)'') value, with negligible influence from 
other matrix elements, including (0)^||iil2||2)''). This is 
because the population of the 4^ state occurs almost ex¬ 
clusively in the two-step i?2-excitation process involving 
the 2'^ state. The population of the 2^ is known very 
well from I-y{2f —>■0)'’) since any significant feeding (i.e. 
from and 2^) can be accounted for. In ^°^Rn, the 
upper limit of {4'^ —)■ 2)'') is used to calculate an upper 
limit for B{E2\ 4f 2)*"), shown in Table IV. 

In the current experiment, the population of the 2^ 
state in ^°^Rn can be considered to occur exclusively via 
a two-step E2 excitation via the 2^ state. The single- 
step process, directly from the ground state, can be as¬ 
sumed to be negligible due to the combination of the 
large energy difference and the small R(if2;2^ —0)'’) 
relative to the R(i?2;2^ —>■ 2^) extracted from the up¬ 
per limit of the branching ratio of 6.9% [47]. Addition¬ 
ally, Ml excitation is calculated to be more than 100 
times weaker than the corresponding E2 between the 
two 2+ states. No complementary data (such as the life¬ 
time, T^+, E2/M1 mixing ratio, ^(2j —>■ 2]*") or conver¬ 
sion coefficient, a(2^ —)■ 2]^)) are available to constrain 
the (2)''IIMlII2J) matrix element and consequently it is 
currently not possible to extract its value. It can how¬ 
ever be shown that the current data set is insensitive 
to the value of the Ml component, and the determi¬ 
nation of (2]''||if2|j2j) is unaffected. The Ml matrix 
element was coupled to the E2 matrix element using 
|5(2j —>■ 2]^) I = 1.1, by comparison to known values in 
the region. 

All of the data for ^°^Rn, collected with both Ag and 
^^°Sn targets, are fitted using the least-squares search 
code, GOSIA [38, 39], in order to fully investigate all po¬ 
tential couplings to unknown matrix elements [35] . In the 
final fit, many matrix elements were coupled, or fixed to 


reasonable values, when the fit was found to be insensi¬ 
tive to their values. The diagonal E2 matrix elements of 
the 4^ and 2^ were coupled to their transitional coun¬ 
terparts, assuming a constant Qq and K — 0 within the 
rigid rotor model. A particular concern with regards to 
correlations is the (4]''||if2||2j) matrix element, which in¬ 
fluences the populations of both the 4^ and 2^ states. It 
was fixed to 0.005 eb in the final fit, although values up 
to 1.5 eb were tested and shown to influence the final 
result at the few percent level, much less than the statis¬ 
tical uncertainty. For the correlated error calculation, it 
was allowed to vary with limits of ±1.5 eb. 

Once the x^ minimum is found, the uncertainties are 
calculated by GOSIA in a two-stage process. At this point, 
all couplings and fixed matrix elements are freed in or¬ 
der to correctly include the influence of correlations to 
unknown matrix elements. Firstly, the diagonal, or un¬ 
correlated, uncertainties on each matrix element are com¬ 
puted by varying it about the minimum until an increase 
in x^ is achieved, satisfying the la condition [39]. At 
the same time, a multi-dimensional correlation matrix 
is built, which is then used in the second step in or¬ 
der to compute the fully correlated errors on each ma¬ 
trix element. It was shown that the {4^1 \\E2\\2'l^) ma¬ 
trix element is insensitive to changes in other transitional 
matrix elements and only very weakly (<< Icr) depen¬ 
dent on (4]^ IIif2114]*"). This leads to an uncertainty on 
B{E2\ d]*" —>■ 2]'') roughly equivalent to the statistical un¬ 
certainty of I'y{4'^ 2^). For (2]^||if2||2j), however, 

the correlations play a much stronger role and the un¬ 
certainty on B{E2;22 -^2'^) is relatively large (see Ta¬ 
ble IV). 

IV. DISCUSSION 

Under the assumption that the quadrupole charge dis¬ 
tribution is uniform and can describe the nuclear shape, 
the deformation can be deduced from the following sum 
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over B{E2) values [48]: 


J2BiE2;0t ^‘2t)=(^ZeRl\ (/3|) , (1) 

i ^ Z 


where Ze is the nuclear charge and Rq = \.2A^I'^ fm. 
From the limit of the 2^ branching ratio in ^°^Rn [47] 
(see Table IV) it is a reasonable assumption that the E2 
transition strength from the ground state is dominated 
by the first-excited 2+ state in these nuclei. Thus, one 
can limit the sum in Equation I to i = 1. The deduced 
deformations are then = 0.099to'oi4 

(/3|)^/^(^°"‘Rn) = 0.120jig;Qjg, indicating a weak deforma¬ 
tion. Another indication of the ground-state deformation 
can be deduced from isotope-shift measurements [18, 19], 
where mean-square charge radii, {r^) a-, is related to the 
deformation (to first order) in the following way: 

+ , ( 2 ) 

where is the mean-square charge radius of a spher¬ 

ical liquid-drop-like nucleus with mass, A [49]. This 
is calculated using the modified liquid-drop model of 
Ref. [50] and the updated parameter set of Ref. [51]. 
Assuming that (/3|)2i2 = 0.062(5), from the Grodzins- 
Raman rule [48] and Equation 2, and using the tabu¬ 
lated A(r^)^^ 2 i 2 values from Ref. [49], (/3 |)a can be de¬ 
duced. All of the derived values for and 

as a function of mass number are plotted in Figure 12 
for comparison. A good level of consistency between the 
two deformation parameters is observed. Furthermore, 
the values obtained from the isotope shift at the heaviest 
masses remain relatively constant apart from the odd- 
even staggering effect. This might be considered as due to 
dynamical effects about a spherical shape, i.e. vibration, 
whereas increase of (/3|) in the lightest isotopes points 
towards an onset of deformation in the ground state. 

A less model-dependent picture of the quadrupole col¬ 
lectivity is the transitional quadrupole moment, Qi, re¬ 
lated to the experimental matrix elements by the follow¬ 
ing relationship: 


Qt{h ^ //) 


jifWEm) I 16 ^ 

(//0201A0) Y 5 (21 f + 1) ’ 


(3) 


where (//020j/i0) is the Clebsch-Gordan coefficient. The 
values deduced from the current experiment are given on 
the level schemes of Figure 13. Here, we can observe that, 
as a function of increasing spin, Qt remains constant in 
^°^Rn as far as the data extends. This can be an indicator 
that these states form a single rotational band, but the 
current level of uncertainty and number of data is not 
enough to make firm conclusions within such a simple 
picture. The alternative and equally simplistic picture 
of an harmonic vibrator gives the relationship between 


Neutron Number, N 

114 116 118 120 122 124 126 


cu 

Q 



FIG. 12. (Color online) Experimental values deduced 

from the B{E2-, 0]^ —> 2^) values measured in this work (black 
circles, “CLX”) and those from isotope-shift measurements 
and liquid-drop model for both the ground (red down trian¬ 
gles, “IS(g)”) and isomeric (red up triangles, “IS(m)”) states. 
The uncertainties on the latter are dominated by the prop¬ 
agation of the uncertainty in the Grodzins-Raman rule [48], 
which is a systematic contribution. The isotope shift values 
are slightly offset from integer A values for clarity of presen¬ 
tation. 


transition strengths of two-phonon {Nph = 2) and one- 
phonon {Nph = 1) states as: 

B{E2; ^ 2 (4) 

BiE2; -t JNj,h=o) 

While this is consistent with the B{E2;if -A 2^) from 
this experiment, it is at odds with the observation of a 
strong B{E 2;22 ‘^i) value, pointing to the fact that 

these nuclei can not be described as simple vibrators (see 
Table IV). For these reasons, comparisons to state-of- 
the-art nuclear models are required to understand the 
behaviour of these nuclei. 

Beyond mean-field calculations have recently been per¬ 
formed for a range of nuclei in this region [52] , having par¬ 
ticular success in describing the electromagnetic matrix 
elements above and below Z = 82 in the polonium [14] 
and mercury [4, 5] isotopes. In these calculations, self- 
consistent mean-field wave functions are generated within 
the Hartree-Fock (HF) -b BCS framework with a Skyrme 
energy-density functional. These are then projected to 
particle number and angular momentum, before being 
mixed by the generator coordinate method (GCM) to 
give physical states. The pure mean-field wave-functions 
are constrained to axial symmetry. While the parame¬ 
ters of the microscopic Skyrme interaction are fitted to 
large sets of data, the extraction of nuclear observables 
from the projected mean-field states can be considered 
parameter free. This is very advantageous when making 
predictions of behaviour where experimental data is not 
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FIG. 13. Comparison of the results of beyond-mean-field calculations and experimental energy levels (italics; in units of keV) 
and transitional quadrupole moments, \Qt\ (in units of eb). The width of the arrows are proportional to \Qt\. Only states up 
to and 2^ are included for clarity of presentation. 


already present. As can be seen in Figure 13, and has also 
been observed in the polonium [14] and mercury [4, 5] iso¬ 
topes, the absolute values of the energy levels predicted 
by the BMF model [52] appear vastly overestimated, but 
the general pattern is reproduced. The prediction of a Oj 
state close in energy to the 2^ state was not able to be 
tested in this experiment. An observation of this state, 
along with its de-excitation branching ratio, would give 
a further test to the model. What is interesting to note 
is the prediction of a very weak 2^ —0]^ decay branch, 
consistent with the observed data, with no need to in¬ 
voke arguments for a forbidden ANp^ = 2 transition in 
the harmonic-vibrator model. 

In Figure 14, the transitional quadrupole moments are 
compared for a range of nuclei extending to A > 194. 
An increase in collectivity for the lighter radon isotopes 
is predicted by an increasing Qt {2f 0]*") value, consis¬ 
tent with the E{2'^) systematics and isotope shift mea¬ 
surements [18, 19]. A more stringent test of this model 
would come from measurements of non-yrast and inter¬ 
band Qt values, which show more significant deviations 
when approaching mid-shell. The current production 
rates at ISOLDE do not allow Coulomb-excitation ex¬ 
periments to be extended to isotopes lighter than ^°°Rn. 
However, with the higher beam energies afforded by HIE- 
ISOLDE [53], multiple-step Coulomb-excitation experi¬ 
ments will be able to provide a complete set of electro¬ 
magnetic matrix elements for the heavier-mass isotopes, 
which may act as a verification of model predictions at 
lower masses. 


V. CONCLUSIONS 

Coulomb excitation of secondary, post-accelerated ra¬ 
dioactive beams of ^°^Rn and ^°^Rn has been performed 
at the REX-ISOLDE facility in CERN. B{E2; 2+ ^ 0]^) 


Neutron Number, N 

_ 108 110 112 114 116 118 



FIG. 14. (Color online) Experimental Qt values in black com¬ 
pared to those predicted by beyond-mean-field calculations in 
red, connected by lines to guide the eye. The upper limits in 
the data are indicated by the downward pointing arrows. 


values have been extracted in both ^°^Rn and ^°'^Rn and 
limits on Qs{2'^) have been determined. In ^°^Rn, popu¬ 
lation of the 2^ and 4]*" states was observed, allowing the 
extraction of B{E2;4:'^ -^2^) and i?(i?2;2j val¬ 

ues in this nucleus. While the excitation energies of the 
observed states in these Rn isotopes coincide with that 
expected of a simple quadrupole vibrator structure, the 
2 ^ —>• 2f transition strength does not support such an in¬ 
terpretation. The results have been compared to recent 
beyond-mean-field calculations [52]. While the energy 
levels seem to be unreasonably expanded, the relative 
behaviour and absolute transitions strengths shows con- 
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sistency between experiment and the model description. 
A more sensitive test of the nuclear shape would come 
from the spectroscopic quadrupole moment, Qs{2^), but 
the precision from this experiment is not sufficient to 
distinguish between oblate-, prolate- and spherical-like 
charge distributions. Extending B(E2) measurements to 
lighter, more exotic nuclei, where shape-coexistence ef¬ 
fects and ground-state deformations are expected to be 
stronger due to the parabolic behaviour of the intrud¬ 
ing structure, would test the BMF description further. 
Observation of a Oj state is still lacking in the light 
radon isotopes. New experiments at higher beam en¬ 
ergy would increase the probability of populating this 
state, should it exist. A future coupling of Miniball 
with the SPEDE electron detector in Coulomb-excitation 
experiments [54, 55] may allow direct detection of the 
if0(0j —>■ 0]*') decay. This will lead not only to a place¬ 
ment of the Oj state in energy, but also to the deter¬ 
mination of the £’2(0J —?► 2[j^)/£0(0j —?► 0]^) branching 
ratio, key to distinguishing between an intruder and a 


phonon structure. Indeed, few-nucleon transfer reactions 
such as {t,p) and (d,p) also have the ability to populate 
such excited 0 + states and could be utilised to elucidate 
their nucleon configuration. 

ACKNOWLEDGMENTS 

We acknowledge the support of the ISOLDE Collab¬ 
oration and technical teams. This work was supported 
by GOA/2010/10 (EOF KULeuven), by the lAP Bel¬ 
gian Science Policy (BriX network P6/23 and P7/12), 
by the U.K. Science and Technology Facilities Coun¬ 
cil (STFC), by the German BMBF under contract Nos. 
05P12PKFNE, 06DA9036I and 05P12RDGIA, by the 
Academy of Finland (Contract No. 131665), by the Eu¬ 
ropean Commission through the Marie Curie Actions 
call PIEFGA-2008-219175 (J.P.). L.P.G. acknowledges 
FWO-Vlaanderen (Belgium) as an FWO Pegasus Marie 
Curie Fellow. 


[1] K. Heyde and J. L. Wood, Rev. Mod. Phys. 83, 1467 

( 2011 ). 

[2] J. Bonn, G. Huber, H.-J. Kluge, L. Kugler, and E. Otten, 
Phys. Lett. B 38, 308 (1972). 

[3] S. Frauendorf and V. Pashkkevich, Phys. Lett. B 55, 365 
(1975). 

[4] N. Bree et ah, Phys. Rev. Lett. 112, 162701 (2014). 

[5] K. Wrzosek-Lipska et ah, Phys. Rev. C (in preparation) 
(2015). 

[6] H. De Witte et ah, Phys. Rev. Lett. 98, 112502 (2007). 

[7] M. D. Seliverstov et ah, Eur. Phys. J. A 41, 315 (2009). 

[8] A. N. Andreyev et ah. Nature 405, 430 (2000). 

[9] P. Rahkila et ah, Phys. Rev. C 82, 011303 (2010). 

[10] T. E. Cocolios et ah, Phys. Rev. Lett. 106, 052503 

( 2011 ). 

[11] M. Seliverstov et ah, Phys. Lett. B 719, 362 (2013). 

[12] R. Julin, K. Helariutta, and M. Muikku, J. Phys. G Nucl. 
Part. Phys. 27, R109 (2001). 

[13] A. Oros et al., Nucl. Phys. A 645, 107 (1999). 

[14] N. Kesteloot et al., Phys. Rev. G (submitted) (2015). 

[15] G. D. Dracoulis, Phys. Rev. C 49, 3324 (1994). 

[16] R. Taylor et al., Phys. Rev. G 54, 2926 (1996). 

[17] R. Taylor et al., Phys. Rev. G 59, 673 (1999). 

[18] W. Borchers et ah, Hyperfine Interact. 34, 25 (1987). 

[19] U. Georg et ah. Laser spectroscopy of ^°°Rn and ^°^Rn 
using alpha detection, in Int. Conf. Exot. Nucl. At. 
Masses ENAM ’95, pages 133-134, 1995. 

U. Georg, Thesis, PhD, Johannes Gutenberg - Univer- 
sitat Mainz, 1997. 

[20] J. Ressler et al., Phys. Rev. C 69, 034317 (2004). 

[21] T. Grahn et ah, EPJ Web Gonf. 63, 01009 (2013). 

[22] D. Dobson et ah, Phys. Rev. G 66, 064321 (2002). 

[23] P. E. Garrett, K. L. Green, and J. L. Wood, Phys. Rev. 
C 78, 044307 (2008). 

[24] J. L. Wood, J. Phys. Conf. Ser. 403, 012011 (2012). 

[25] E. Clement et ah, Phys. Rev. C 75, 054313 (2007). 

[26] L. P. Gaffney et ah. Nature 497, 199 (2013). 


[27] D. G. Jenkins, Nat. Phys. 10, 909 (2014). 

[28] L. Penescu, R. Catherall, J. Lettry, and T. Stora, Rev. 
Sci. Instrum. 81 , 02A906 (2010). 

[29] B. H. Wolf et al., Nucl. Instrum. Meth. B 204, 428 
(2003). 

[30] F. Wenander, J. Instrum. 5, C10004 (2010). 

[31] D. Voulot et ah, Nucl. Instrum. Meth. B 266, 4103 
(2008). 

[32] N. Warr et al., Eur. Phys. J. A 49, 40 (2013). 

[33] S. Zhu and F. Kondev, Nucl. Data Sheets 109 , 699 
(2008). 

[34] J. F. Ziegler, M. D. Ziegler, and J. P. Biersack, Nucl. 
Instrum. Meth. B 268, 1818 (2010). 

[35] M. Zieliriksa, L. P. Gaffney, and K. Wrzosek-Lipska, (in 
preparation) (2015). 

[36] N. Bree et al., Nucl. Instrum. Meth. B (2015). 

[37] O. Helene, Nucl. Instrum. Methods 212, 319 (1983). 

[38] T. Czosnyka, D. Cline, and C. Y. Wu, Bull. Am. Phys. 
Soc. 28 (1983). 

[39] D. Cline et ah, Gosia User Manual for Simulation 
and Analysis of Coulomb Excitation Experiments, Gosia 
Steering Committee, 2012. 

[40] T. Kibedi, T. Burrows, M. Trzhaskovskaya, P. Davidson, 
and C. Nestor, Nucl. Instrum. Meth. A 589, 202 (2008). 

[41] R. Robinson, P. McGowan, P. Stelson, and W. Milner, 
Nucl. Phys. A 150, 225 (1970). 

[42] M. Zielihska et al., HIL Annual Report; Coulomb exci¬ 
tation ^°^Ag, Technical report. Heavy Ion Laboratory, 
University of Warsaw, Warsaw, Poland, 2009. 

[43] M. Throop, I. Hall, I. Naqib, D. Thomas, and B. Wake¬ 
field, Phys. Lett. B 41, 585 (1972). 

[44] J. Blachot, Nucl. Data Sheets 107 , 355 (2006). 

[45] T. Miller and M. Takeda, Nucl. Phys. A 221 , 392 (1974). 

[46] M. Loiselet, O. Naviliat-Cuncic, and J. Vervier, Nucl. 
Phys. A 496, 559 (1989). 

[47] S. Gillespie, private communications (2014). 



12 


[48] S. Raman, C. Nestor, and P. Tikkanen, At. Data Nucl. 
Data Tables 78, 1 (2001). 

[49] E. W. Otten, Treatise Heavy-Ion Sci. 8, 517 (1989). 

[50] W. D. Myers and K.-H. Schmidt, Nucl. Phys. A 410, 61 
(1983). 

[51] D. Berdichevsky and F. Tondeur, Zeitschrift fiir Phys. A 
Atoms Nucl. 322, 141 (1985). 


[52] J. M. Yao, M. Bender, and P.-H. Heenen, Phys. Rev. C 
87, 034322 (2013). 

[53] M. Lindroos, P. Butler, M. Huyse, and K. Riisager, Nucl. 
Instrum. Meth. B 266, 4687 (2008). 

[54] J. Konki et ah, EPJ Web Conf. 63, 01019 (2013). 

[55] J. Pakarinen et ah. In-beam electron spectroscopy at 
HIE-ISOLDE, Technical Report CERN-INTC-2010-040. 
INTC-I-107, CERN, Geneva, 2010. 


